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OBSERVATION  OF  BREWSTER  ANGLE  LIGHT  SCATTERING  FROM  AIR 

BUBBLES  RISING  IN  WATER 
Abstract 


by  Stefan  M.  BSumer,  M.S. 
Washington  State  University 
August  1988 


Chair:  Philip  L.  Marston 

Air  bubbles  rising  freely  through  distilled  water  were  illuminated  by  laser  light  with 
incident  polarization  of  the  electric  field  parallel  to  the  plane  of  incidence.  The  radii  a  of  the 
bubbles  investigated  were  in  the  range  50  |im  <  a  <  100  (am.  At  Brewster's  angle  of 
incidence  is  =  arctan  (m)  =  36.9°,  the  Fresnel  reflection  coefficient  vanishes,  where  m  is 
the  relative  index  of  refraction  for  a  water  -  air  interface,  m  =  0.75.  Therefore  a  minimum  in 
scattered  intensity  should  be  observed  at  the  Brewster  scattering  angle 
0b  =  180°  -  2ie  =  106.2°.  Furthermore  interference  patterns  caused  by  different  classes  of 
rays  should  change  drastically  for  angles  around  Brewster  scattering  angle  0b  Experiments 
were  done  to  detect  these  two  phenomena  by  photographing  the  scattered  light.  The 
scattered  intensity  as  a  function  of  angle  was  inferred  by  scanning  the  photographs  with  a 
microdensitometer.  In  order  to  confirm  the  experimental  results,  they  were  compared  with 
Mie  theory  calculation  for  scattering  from  spherical  bubbles.  Features  of  the  scattering 
pattern  due  to  ray  interference  could  be  understood  using  geometrical  optics.  The 
experiment  confirmed  that  the  observed  features  of  the  scattering  pattern  near  0b  are  those 
predicted  of  a  bubble  having  a  clean  surface.  Since  it  has  been  previously  calculated  that  the 
scattering  pattern  would  be  significantly  different  near  0b  if  the  surface  of  the  bubble  is 
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coated  by  contamination,  it  is  plausible  that  such  a  coating  could  be  characterized  optically. 
It  is  thought  that  the  acoustics  and  dynamics  of  microbubbles  in  the  ocean,  or  in  other 
environments,  may  be  affected  by  such  coatings.  The  present  observations  could  be 
compared  with  observations  of  the  scattering  from  coated  bubble"  if  such  observations 
become  available. 
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Chapter  1 

Introducflon 

Light  scattering  from  small  particles  has  been  of  interest  for  a  long  time  for  various 
reasons.  One  particular  class  of  scatterers  has  attracted  some  recent  attention:  air  bubbles  in 
dielectric  media.  Some  research  has  been  done  in  this  field  experimentally  and 
theoretically1'11.  However  there  are  still  areas  in  which  little  or  no  research  was  done.  One 
of  these  areas  is  Brewster  angle  light  scattering.  Investigations  on  this  particular  topic  in 
scattering  will  be  presented  in  the  following  work.  Experimental  data  will  be  shown  and 
compared  with  theoretical  results.  The  motivation  for  this  work  is  twofold:  interest  in  the 
phenomena  itself  and  possible  applications  for  identifying  and  classifying  microbubbles  in 
water7. 

First  the  theory  of  Brewster's  angle  will  be  described  and  a  model  will  be  presented 
which  was  used  in  understanding  the  phenomena  and  in  calculating  theoretical  data. 
Chapter  three  describes  the  experiment,  and  in  chapters  four  and  five  results  will  be 
presented  and  discussed. 


Chapter  2 

Theory 


If  an  electromagnetic  light  wave  propagates  through  a  medium  and  hits  an  interface  with 
another  medium,  two  phenomena  occur:  refraction  and  reflection.  Refraction  and  reflection 
are  qualitatively  described  by  Snell's  law  and  the  law  of  reflection.  The  direction  of  the 


refracted  light  is  given  by  Snell's  law 
ni  sin  ii  =  nj'  sin  ij' 


(1) 


2 


where  ni  and  ni'  are  the  indices  of  refraction  for  media  1  and  2,  ii  and  ij'  are  angle  of 
incidence  and  angle  of  refraction. 

The  direction  of  the  reflected  light  is  determined  by  the  law  of  reflection 
ii  =  if 

where  ir  is  the  angle  of  reflection. 

In  addition  there  is  the  boundary  condition  that  both  the  reflected  and  refracted  ray  are  in  the 
same  plane  as  the  incident  ray. 

If  we  are  not  only  interested  in  the  direction  of  the  light  but  also  in  the  magnitude  of  the 
reflected  and  refracted  light  wave,  we  have  to  use  Fresnel's  reflection  and  transmission 
coefficients. 

Eout  =  r  Ein 

The  outgoing  electric  field  amplitude  Eom  equals  the  Fresnel  coefficient  r  times  the  incident 
electric  field  amplitude  £„.  The  Fresnel  coefficients  depend  on  ii  and  ii  and  also  on  the 
incident  polarization  of  the  light.  Therefore  we  have  four  such  coefficients,  two  for 
reflection  and  two  for  transmission.  The  following  list  is  taken  from  Klein-Furtak19. 


tan(ii  -  ii') 
p  tan(ii  +  i  i ') 


Fresnel  reflection  coefficient,  electric  field  parallel  polarized 


■  sin(ii-ii') 
*  ’  sin(ii  +  ii ') 


Fresnel  reflection  coefficient,  electric  field  perpendicular  polarized 


t  -  2  sin  ii  cos  ii  Fresnel  transmission  coefficient 

P  sin(ii  +  u  )  cos(m  -  i i ' ) 

electric  field  parallel  polarized 


2  sin  ij’  cos  ii 
ts  =  sin(i,+  l,1) 


\ 


Fresnel  transmission  coefficient, 
electric  field  perpendicular  polarized 


3 


Parallel  and  perpendicular  always  refer  to  the  plane  of  incidence,  which  in  this  case  is 
identical  with  the  plane  of  the  paper  in  the  following  figures. 

As  we  examine  the  reflection  coefficient  for  parallel  polarized  light,  we  notice  that  if 
ii  +  ii'  =  90°  (2) 

Mmr  i;  tan(ii-ii')  „ 

limrp-Um^r— j-r-ry  _0 

ii+ii'_»90° 

The  reflection  coefficient  for  parallel  polarized  light  incident  rp  goes  to  zero. 

Applying  this  special  condition  (2)  to  Snell's  law  (1)  we  find 
ni  sin  ii  =  m'  sin  ii*  =  ni'  sin(90°  -  ii)  =  ni'  cos  ii 


tan  it  -jjJ-  (3) 

So  condition  (2)  leads  to  the  definition  of  a  special  angle  of  incidence.The  angle  defined 
through  equation  (3)  is  called  Brewster's  angle  is,  named  after  Sir  David  Brewster  (1781  - 
1868),  who  derived  this  law  in  1815  empirically21.  At  Brewster's  angle  no  parallel 
polarized  light  is  reflected,  all  parallel  polarized  light  is  transmitted. 


Now  we  want  to  apply  this  general  relation  to  the  specific  case  of  the  water  -  air 
interface  of  an  air  bubble  freely  rising  through  water.  First  of  all  the  size  range  of  interest 
for  bubbles  in  this  experiment  is  50  (im  <  a  <  100  |im,  where  a  is  the  bubble  radius.  As 
shown  in  previous  work4-11-13  microbubbles  of  that  size  can  be  well  modeled  as  spheres. 
The  sphere  model  will  be  used  throughout  all  calculations.  Since  the  radii  of  the  bubbles  of 
interest  are  a  lot  larger  than  the  wavelength  of  the  light,  we  can  treat  all  interfaces  of  the 
bubbles  locally  as  plane  interfaces,  and  apply  the  above  derived  formulas  directly. 

If  a  light  wave  is  incident  on  the  water  -  air  interface,  as  shown  in  figure  1,  it  makes  a 
local  angle  of  incidence  i  with  the  local  normal  defined  by  the  extension  of  a  line  from  the 
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Figure  1:  Geometry  of  light  rays  incident  onto  an  air  bubble  in  water 
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center  of  the  bubble  to  the  point  of  incidence.  Applying  Snell's  law  and  the  law  of 
reflection  we  also  mark  the  angle  of  reflection  4  which  has  the  same  magnitude  as  i,  and  the 
angle  of  refraction  i'.  The  scattering  angle  0s  is  defined  by  the  forward  scattering  direction 
as  indicated  in  figure  1.  Lastly  there  is  the  impact  parameter  s,  which  is  defined  by  the 
distance  of  the  incident  light  ray  to  the  optical  axis,  where  the  optical  axis  goes  through  the 
center  of  the  bubble.  Now  that  we  have  named  all  components  we  need  to  completely 
describe  the  geometry  of  our  problem.  The  next  step  is  to  calculate  Brewster’s  angle  in  case 
of  water  to  air  transmission.  Using  equation  (3)  we  calculate 

iB  =  arc  tan  m  =  arc  tan  0.75  =  36.9°  (4) 

which  is  Brewster's  angle  for  light  going  from  water  into  air.  In  this  equation  m  is  the 
relative  index  of  refraction  defined  by 

«> 

As  the  numerical  value  m  =  0.75  was  taken.  Although  the  ratio  na/nw  of  the  actual  values 
for  na  and  nw  at  514.5  nm  is  slightly  different  (na= 1.0002929,  nw=  1.334, 
m=0. 74984) 15,12,  m=0.75  has  proven  to  be  a  very  good  approximation  in  all 
calculations1,4.  Going  back  to  figure  (1)  we  find  the  scattering  angle  to  be 

0s  =  180° -  2 i  (6) 

Using  the  Brewster  angle  incident  (equation  4),  we  can  compute  the  Brewster  scattering 
angle  0B. 

0B  =  180°  -  2  iB  =  106.2° 

Besides  understanding  how  to  model  the  reflection  of  light  of  the  surface  of  a  bubble 
we  also  have  to  know  what  happens  to  the  light  rays  that  are  transmitted  into  the  air  bubble. 
It  is  of  particular  interest  whether  light  with  an  incident  angle  other  than  Brewster's  angle  is 
transmitted  out  of  the  bubble  at  the  Brewster  scattering  angle  0b  after  internal  reflections . 
This  is  because  if  light  is  emitted  near  the  same  direction  as  Brewster  scattering  angle,  we 
expect  to  see  modulated  interference  patterns,  because  the  observed  light  comes  from  the 


same  source  but  travels  paths  of  different  lengths.  The  interference  patterns  then  pro  vide 
information  about  the  bubble.  To  gain  information  about  all  light  rays  scattered  in  the  same 
direction  we  apply  the  following  generalized  formula  for  finding  scattering  angles  from 
bubbles11-17: 

0s  =  h(2pir-2i- Jt(p-2g-l))  (7) 

0s,  i,  ir  are  as  before  and 

p  number  of  internal  chords  for  the  light  ray 
g  number  of  trips  around  the  bubble  for  internally  reflected  light  rays 

h  +1  for  rays  entering  above  the  optical  axis 

-1  for  rays  entering  below  the  optical  axis 

For  the  special  case  of  p  =  0,  g  =  0  we  get  the  result  of  equation  (6),  Brewster's  scattering 
angle.  If  we  scan  through  different  incident  angles  from  0°  to  the  critical  angle  ic=aicsin(m), 
we  see  whether  0s  goes  through  Brewster's  scattering  angle.  We  do  not  get  a  contribution 
from  one-chord  rays,  but  one  two-chord  ray  and  two  three-chord  rays  contribute  (table  1). 
However  finding  the  exact  angle  of  incidence  for  those  rays  is  a  non  trivial  problem,  since 
Snell's  law  (1)  and  equation  (7)  have  to  be  solved  simultaneously.  That  leads  to  a 
transcendental  equation  which  we  solved  numerically  by  rewriting  equations  (4)  and  (6)  as 
0b  =  2  arc  tan  (nr1)  (8) 

and  finding  the  roots  of  f(i)  =  0,  where 

f(i)  *  2  arctanlnr1)  -  h(2pir  -  2i  -  k  (p  -2g  -1))  (9) 

The  result  is  shown  in  table  2. 

With  this  information  we  are  now  able  to  calculate  the  impact  parameters  s  for  the  relevant 
rays 

s  =  asini  (10) 

Knowing  the  impact  parameters  and  applying  Snell's  law  enables  us  to  draw  a  ray  diagram 
for  Brewster  angle  scattering  (figure  2).  In  the  notation  used,  the  first  number  stands  for 


Table  1:  Numerical  search  for  Brewster  scattering  angle  for  different  types  of  rays 

Number  of  chords  in  the  bubble 
2 

Number  of  trips  around  the  bubble 
0 


gle  of  incidence 

Scattering  angl 

1 

176.6665 

2 

173.3316 

3 

169.9943 

4 

166.6531 

5 

163 . 3069 

6 

159.9541 

7 

156.5936 

8 

153.2239 

9 

149.8436 

10 

146.4511 

11 

143 . 0450 

12 

139.6236 

13 

136.1853 

14 

132.7281 

15 

129.2503 

16 

125.7499 

17 

122.2246 

18 

118.6722 

19 

115.0902 

20 

111.4759 

21 

107.8265 

22 

104.1388 

23 

100.4093 

24 

96.6344 

25 

92 . 8097 

26 

88.9307 

27 

84.9923 

28 

80.9887 

29 

76.9133 

30 

72 .7587 

31 

68 . 5167 

32 

64 . 1775 

33 

59 .7298 

34 

55. 1606 

35 

50.4544 

36 

45.5925 

37 

40.5527 

39 

35.3073 

39 

29.8213 

40 

24.0521 

4  1 

17 . 9406 

42 

11. 4091 

4  3 

4.3477 

*4 

-3.4063 

i  D 

-12.1151 

46 

-22.2425 

«  n 

/ 

-34.7335 

43 

-52 . 9396 

Number  of  chords  in  the  bubble 
3 

Number  of  trips  around  the  bubble 
0 


of  incidence 

Scattering  angle 

1 

353.9997 

2 

347.9975 

3 

341.9915 

4 

335.9797 

5 

329.9603 

6 

323.9312 

7 

317.8904 

8 

311.8358 

9 

305.7654 

10 

299.6767 

11 

293.5675 

12 

287.4355 

13 

281.2779 

14 

275.0922 

15 

268.8755 

16 

262.6248 

17 

256.3369 

18 

250.0083 

19 

243.6353 

20 

237.2139 

21 

230.7397 

22 

224.2082 

23 

217 . 6140 

24 

210.9516 

25 

204 .2146 

26 

197.3961 

27 

190.4885 

28 

183.4830 

29 

176 . 3699 

30 

169 . 1381 

31 

161.7751 

32 

154 .2662 

33 

146 . 5947 

34 

138.7409 

35 

130.6815 

36 

122 . 3888 

37 

113.8291  ! 

38 

104.9609  ! 

39 

95.7327 

40 

86.0782 

4  1 

75.9109 

42 

65.1136 

4  3 

53.5215 

44 

40.8898 

45 

26.8273 

46 

10.6362 

47 

-9.1828 

48 

-37.4844 

WK0 


number  of  chords  in  the  bubble,  the  second  number  tells  how  many  times  a  ray  crosses  the 
optical  axis. 
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Since  we  found  several  rays  contributing  to  the  scattering,  we  expect  to  see  interference 
patterns  in  the  Brewster  angle  region.  To  characterize  the  interference  we  use  methods  of 
geometrical  optics  and  calculate  the  separation  of  adjacent  extreme  A016. 


A0  =  ^f  (11) 

where  Xw  is  the  wavelength  of  the  incident  light  in  water  and  d  the  distance  at  the  bubble 
exit  plane  between  two  interfering  rays.  Looking  at  figure  (2)  we  find  by  symmetry  that  the 
separation  d  of  principal  rays  in  the  exit  plane  is  given  by  the  sum  of  their  impact 


parameters  so  and  Si 

d  =  so+  Sj 

Therefore  we  get  the  final  result  as1-4 


So+S]  v  ' 

Note  that  the  greater  d  is,  the  smaller  A0  becomes.  This  means  high  frequency  oscillations 
are  caused  by  widely  spaced  rays. 

The  complete  general  solution  for  a  plane  wave  scattered  of  a  dielectric  sphere  was  given 
by  G.  Mie  in  1908  20.  But  only  recently  with  the  advent  of  computers  it  was  possible  to 
calculate  fairly  easily  scattered  intensities  over  a  range  of  angles  for  bubble  sizes  much 
larger  than  the  wavelength  of  the  used  light .  Figure  (3)  shows  the  Mie  computation  result 
for  a  bubble  of  typical  size  in  this  experiment.  The  bubble  is  characterized  by  the  size 
parameter  ka,  where  a  is  the  bubble  radius  and  k=2  rr.9.w  .  All  uuuiMties  in  the  Mie 


calculations  are  normalized  with  respect  to  a  perfectly  reflecting  sphere  of  the  same  size. 

For  the  computations  the  original  MIEVO  code  by  W.J.  Wiscombe18  modified  by 
S.C.Billette9  and  C.E.Dean7  was  taken.  We  clearly  see  the  expected  interference  pattern 
and  also  a  broad  local  minimum  in  intensity.  A  similar  result  is  expected  to  be  obtained  in 
the  experiment. 

To  describe  the  origin  of  the  scattering  pattern  around  Brewster's  angle  caused  by 
interference,  we  now  calculate  d/a  for  Brewster  angle  scattering  using  equations  10,12. 
Since  this  quantity  is  proportional  to  A0  we  are  able  to  get  first  information  on  what 
frequencies  to  expect.  The  result  of  d/a  calculated  for  all  theoretical  possible  interferences  is 
shown  in  table  3.  We  notice  that  d/a  is  approximately  the  same  for  (0,0),(2,1)  and 
(3,1), (3,2)  interference.  Therefore  we  do  not  expect  to  see  a  separate  (3,1),(3,2) 
interference.  A  similar  situation  occurs  between  the  (0,0), (3,2)  and  (2,1),(3,1) 
interferences.  Hence  a  distinct  (2,1), (3,1)  frequency  might  also  be  hard  to  observe. 

Besides  this  observation  we  also  have  to  consider  the  amplitude  of  the  contributing  rays 
in  the  far  field  in  order  to  make  quantitative  predictions  about  which  interference  patterns 
we  expect  to  see.  Concerning  the  amplitude  in  the  far  field,  several  considerations  have  to 
^  "  ide  ” .  In  our  case  we  choose  to  pay  attention  only  to  one  very  important  fact:  how 
interfaces  each  light  ray  hits.  Every  time  light  hits  an  interface,  the  intensity  of  the 
remaining  light  will  be  reduced  because  both,  reflection  and  transmission  occur.  Therefore 
(0,0)  and  (2,1)  are  ordinarily  dominant  rays.  Since  they  are  widely  spaced  they  cause  the 
easy  to  detect  high  frequency  (closely  spaced)  oscillations1-4.  Around  106°  the  (0,0)  ray 
goes  through  a  local  minimum  in  intensity.  At  Brewster’s  angle  all  parallel  polarized  light  is 
transmitted  into  the  bubble,  none  reflected.  Because  of  that  the  scattered  intensity  also  goes 
through  a  local  minimum  and  the  high  frequency  oscillation  vanishes  and  then  reappears. 
The  remaining  structure  around  106°  is  caused  by  other  interferences  than  those  involving 
the  (0,0)  ray.  The  fringes  around  106°  are  due  to  (2,1), (3,1)  or  (2,1), (3, 2)  interferences. 
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interference  between 

d/a  =  sin  i  +  sin  j 

(0,0) 

(2,1) 

0.965300146 

(0,0) 

(3,1) 

1.213676858 

(0,0) 

(3,2) 

0.300773020 

(2,1) 

(3,1) 

0.248376712 

(2,1) 

(3,2) 

0.664527126 

(3,1) 

(3,2) 

0.912903839 

Table  3;  d/a  for  all  possible  interferences  with  p  <  3 
i  angle  of  incidence  for  the  first  ray 
j  angle  of  incidence  for  the  second  ray 
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Furthermore  we  notice  a  lower  frequency  modulation  outside  the  Brewster  angle  scattering 
region.  Because  of  the  low  frequency  it  is  most  likely  caused  by  (0,0),(3,2)  interference, 
since  these  rays  are  closely  spaced. 

The  experimental  goal  is  to  find  some  of  these  effects  for  air  bubbles  rising  freely 
through  clean  water 

i)  the  disappearance  of  the  fringe  pattern 

ii)  a  broad  local  minimum  in  scattered  intensity 

iii)  expected  interference  patterns. 


Chapter  3 


ixoenment 


The  basic  set  up  which  was  used  throughout  the  experiment  is  shown  in  figure  4.  The 
light  was  provided  by  a  Spectra  Physics  argon  ion  laser  operating  at  5 14.5  nm  in  air  with  a 
gaussian  beam  of  1.25  mm  diameter  at  the  1/e2  point  A  high  quality  ellipsometer  grade 
polarizer  was  used  to  get  a  very  accurate  polarization  because  Brewster's  angle  occurs  only 
for  parallel  polarized  light  as  shown  earlier.  The  next  component  in  the  light  path  was  an 
aperture.  It  served  to  clean  the  laser  beam,  i.e.  blocking  out  reflections  from  previous 
devices  in  the  beam  and  eliminating  low  frequency  modulations.  Behind  the  aperture  is  the 
scattering  chamber.  It  was  specially  built  for  this  experiment  and  has  some  specifications 
that  arc  worth  mentioning.  The  following  description  of  the  chamber  refers  to  figures  5  -  6. 

The  scattering  chamber  is  built  out  of  aluminum.  It  has  three  windows,  an  entrance 
window,  an  exit  window,  and  an  observation  window.  Through  the  entrance  window  at 
the  front  side  of  the  chamber  the  laser  beam  enters  the  chamber.  It  is  oriented  so  that  the 
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Ray  diagram 


for  Brewster  angle  light  scattering 


window 


incident  light  is  perpendicular  to  the  window  in  order  to  minimize  refraction.  The  exit 
window  at  the  back  of  the  chamber  is  tilted  2°  with  respect  to  the  entrance  window.  The 
reason  for  this  is  to  avoid  illumination  of  the  bubble  by  the  reflection  from  the  rear 
window.  Consider  the  case  of  parallel  entrance  and  exit  windows.  We  can  calculate  the 
scattering  angle  of  the  reflected  light  with  the  bubble.  This  gives  the  following  result  for  the 
new  scattering  angle  0s 

es'  =  180°  -  0B  =  1 80°  -  106°  =  74° 

This  new  angle  0s'  is  very  close  to  the  region  of  critical  angle  scattering3.  Although  the 
laser  beam  is  only  partially  reflected  from  the  exit  window,  the  effect  of  critical  angle 
scattering  from  the  reflected  beam  is  still  a  lot  stronger  than  Brewster  angle  scattering  from 
the  original  beam  in  that  region.  So  all  one  would  see,  provided  parallel  entrance  and  exit 
windows,  would  be  the  scattering  pattern  due  to  critical  angle  scattering  which  would 
completely  dominate  the  Brewster  angle  effect.  (A  critical  angle  scattering  pattern  was 
actually  observed  with  a  preliminary  version  of  the  apparatus  in  which  the  exit  window  was 
not  tilted.)  Therefore  the  rear  window  was  tilted  so  that  the  reflected  beam  goes  off  in  a 
different  direction  and  misses  the  bubble.  Finally  we  have  the  observation  window  at  one 
side  of  chamber.  This  window  is  oriented  so  that  it  is  perpendicular  to  the  scattered  light 
coming  of  at  106°,  the  calculated  Brewster  scattering  angle.  All  windows  were  high  quality 
optical  windows.  The  entrance  and  observation  window  had  antireflection  coatings  on  the 
outside.  That  helped  to  reduce  unwanted  reflected  light  off  the  windows.  The  whole  body 
of  the  chamber,  inside  and  outside,  was  anodized  black  in  order  to  minimize  reflections  of 
the  walls.  At  the  bottom  of  the  chamber  a  small  pipe  is  mounted.  At  the  bottom  of  this  pipe 
a  hollow  needle  is  placed  which  is  connected  to  an  air  supply.  By  regulating  the  air  through 
the  needle  bubbles  were  created.  While  rising  through  this  small  pipe  the  bubble  stabilizes, 
since  almost  all  of  the  internal  size  oscillations  damp  out  in  a  few  parts  of  a  second.  On  top 
of  the  chamber  body,  concentric  with  the  bottom  pipe,  there  is  another  pipe,  larger  in 


diameter.  This  pipe  provides  time  to  prepare  for  capturing  and  measuring  the  bubble,  since 
it  takes  a  couple  of  seconds  for  the  bubble  to  rise  to  the  top  of  that  pipe  after  the  scattering 


event. 

The  last  component  used  in  the  experiment  was  a  camera  for  recording  the  scattering 
event  It  was  mounted  in  front  of  the  observation  window.  In  order  to  reduce  background 
scattering  from  the  water  and  from  the  various  interfaces,  a  polarization  filter  was  placed  in 
front  of  the  camera  lens.  The  polarization  filter  was  set  to  pass  light  polarized  parallel  to  the 
scattering  plane. 

Because  the  effect  we  are  looking  for  occurs  only  at  a  specific  angle,  it  is  essential  that 
all  components  remain  in  the  same  position  relative  to  each  other.  To  make  sure  of  that  all 
components  were  clamped  to  a  rigid  reference  frame. 

Doing  the  actual  experiment  consists  of  three  parts: 

i)  angle  calibration 

ii)  background  determination 

iii)  data  taking 

3.2  Angle  calibration 

The  first  step  in  each  run  of  the  experiment  is  angle  calibration,  since  the  effect  we  are 
looking  for  occurs  only  at  a  specific  angle.Therefore  we  need  to  identify  angles  precisely 
on  the  recorded  scattering  pattern.  In  order  to  get  such  a  precise  calibration,  we  had  to  add  a 
couple  of  components  to  the  basic  set  up  (figure  7).  First  of  all  we  need  a  reference 
direction.  For  this  the  back  scattering  direction  was  taken.  To  find  the  exact  back  scattering 
direction  a  beamsplitter  was  inserted  in  the  light  path.  Then  a  retroreflector  was  placed 

|  behind  the  beamsplitter.  The  beam  coming  out  of  the  retroreflector  is  reflected  at  45°  off  the 

\ 
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1  beamsplitter  and  then  focused  to  a  point.  This  point  now  defines  the  backwards  direction. 
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Scattering  chamber 


retroreflector 


In  the  next  step  a  thin  glass  slide,  serving  as  a  simple  beam  splitter,  is  placed  in  the 
scattering  chamber  (figure  8).  The  glass  slide  is  mounted  to  the  end  of  a  rod  which  is 


connected  to  a  goniometer.  By  rotating  the  glass  slide  to  known  amounts  indicated  by  the 
goniometer,  the  incident  light  beam  can  be  reflected  in  all  wanted  directions.  But  before 
applying  this  technique  for  angle  calibration,  the  goniometer  had  to  be  calibrated.  We  rotate 
the  glass  slide  so  that  the  reflected  ray  of  the  glass  slide  is  focused  into  the  point  defined 
earlier  by  the  retroreflector.  The  reading  the  goniometer  shows,  now  defines  the  reference 
backwards  direction  on  the  goniometer. 

Finally  the  calibrated  goniometer  was  used  to  mark  reference  points  of  exactly  known 
angle  on  photographic  film.  The  106°  scattering  angle  direction  was  dialed  on  the 
goniometer.  Several  dots  around  106°  are  marked  on  one  multiple  exposed  frame  of  the 
film.  For  each  dot,  the  exact  angle  is  known.  Typically  about  six  calibration  points  were 
taken.  With  this  information  we  were  able  to  tell  the  scattering  angle  of  every  point  on  the 
data  negatives.  This  was  done  by  using  a  linear  extrapolation  to  determine  the  angle  for 
each  point  on  the  picture.  One  edge  of  the  negative  always  served  as  initial  reference  for 
applying  the  angle  extrapolation  to  different  negatives. 

The  whole  angle  calibration  procedure  is  similar  to  that  introduced  by  Langley  and 
Marston4.  It  also  automatically  accounts  for  the  effects  of  refraction  at  the  water  -  window  - 
air  interface. 

3.3  Background  determination 

Besides  the  desired  scattering  from  the  bubble,  mere  is  always  unwanted  background 
scattering  which  is  also  recorded  on  the  film.  In  order  to  compare  the  data  with  theoretical 
results,  we  have  to  determine  the  background  and  subtract  it  out  of  the  recorded  data.  A 
picture  without  bubble  in  the  beam  was  taken.  This  represents  all  the  background  scattering 
that  is  present.  Since  the  background  is  not  a  function  of  time,  one  background  picture  on 


each  roll  of  film  can  represent  the  background  for  all  data  in  that  particular  run  of  the 
experiment.  Such  a  negative  was  analyzed  and  a  background  function  was  found  as 
discussed  in  detail  in  a  later  section. 


3.4  Data  taking 

The  actual  process  of  taking  data  is  a  fairly  simple  one.  An  air  bubble  is  created  with  the 
needle  at  the  bottom  of  the  lower  pipe.  Then  the  bubble  is  watched  rising  freely  through  the 
water.  When  it  is  about  to  enter  the  light  beam,  the  camera  is  triggered  by  hand.  The  shutter 
was  closed  as  soon  as  the  bubble  left  the  beam  again.  The  exposure  time  was  measured 
with  a  simple  clock  triggered  through  the  camera.  After  that,  the  bubble  was  captured  under 
a  glass  slide,  which  was  put  on  top  of  the  upper  pipe.  While  sitting  under  the  glass  slide  the 
diameter  of  each  bubble  was  measured  using  a  measuring  microscope. 

During  the  whole  process  the  camera  was  focused  on  infinity  so  that  the  Fraunhofer  far 
field  approximation  yields.  In  different  experiments  a  Nikon  85  mm  lens  and  a  Nikon  135 
mm  lens  together  with  a  Nikon  F2  camera  were  used.  For  photographic  film,  Kodak  Tri  X 
pan  400  ASA  was  used.  The  laser  delivered  typically  300  mW  power  output.  Throughout 
the  experiment  distiled  water  was  used  to  reduce  the  amount  of  impurities  which  would 
otherwise  tend  to  coat  the  bubble7-9.  Furthermore  distilled  water  helps  reducing 
background  scattering  from  small  particles  as  well  as  getting  as  close  as  possible  to  the 
assumed  refractive  index12. 


: 


26 


Chapter  4 

Results 


4.1  Photographic  data 

In  figure  9  a  typical  photograph  of  a  scattering  event  is  shown.  We  are  able  to  recognize 
the  high  frequency  oscillation  in  intensity  resulting  from  the  predicted  interference  patterns. 
Furthermore  a  disappearing  of  the  fine  structure  in  the  middle  of  the  picture  can  be  noticed. 
This  is  first  experimental  evidence  of  Brewster  angle  scattering.  But  in  order  to  get  more 
detailed  information  from  the  picture  concerning  the  spacing  of  the  fringes,  relative 
intensities  and  angles,  the  data  from  the  photographic  negatives  have  to  be  processed. 

4.2  Data  processing 

After  developing  the  film  a  digital  intensity  profile  was  taken  with  a  microdensitometer. 
To  get  good  results,  the  following  procedure  was  applied.  Five  scanlines  in  horizontal 
direction  in  the  middle  of  the  negative  were  taken.  The  scanlines  were  0.2  mm  apart  from 
each  other.  The  scanning  was  done  with  the  scanning  optics  defocused.  The  reason  for  that 
is  that  the  microdensitometer  resolves  structures  at  about  the  grain  size  of  the  negative.  If  a 
negative  is  scanned  somewhat  out  of  focus,  optical  averaging  in  the  horizontal  direction  is 
done.  This  is  done  to  reduce  grain  noise  and  emphasize  the  real  structure  of  the  data.  Even 
when  defocused,  the  resolution  of  the  microdensitometer  greatly  exceeds  the  spacing  of  the 
fine  fringes.  Other  sources  of  random  noise  are  imperfections  of  negative  and  scanning 
optics  and  dust.  These  effects  were  reduced  by  taking  the  five  scanlines  and  by  forming 
mathematical  average.  Usually  760  data  points  were  taken  across  the  negative.  Considering 
around  600  data  points  in  a  relevant  angle  range  of  about  12°,  we  end  up  with  a  stepsize  ot 
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around  0.C2  degrees/data  point.  The  exact  stepsize  depends  on  the  camera  lens  which  was 
used  to  record  the  specific  negative. 

Besides  reducing  the  noise  and  bringing  out  the  structure  of  the  data  as  much  as 
possible,  there  is  another  problem.  In  order  to  compare  experimental  data  with  theoretical 
calculations,  we  ideally  would  like  to  know  absolute  intensities  for  the  recorded  scattering 
pattern.  The  microdensitometer  puts  out  only  a  relative  intensity  profile.  Therefore  we  have 
to  use  an  empirically  obtained  response  curve  of  the  photographic  film.  This  curve  is  called 
the  Hurter  -  Driffield  (H&D)  curve14.  In  this  curve  the  logarithm  of  the  exposure  logE  is 
plotted  versus  the  photographic  density  D.  The  exposure  is  defined  as 

E  =  J  T  (14) 

with  J  as  the  intensity  incident  during  the  exposure  time  T, 

and  the  photographic  density  is  determined  by 

D  =  log(-)  (15) 

x 

where  I  =  with  Ittans  being  the  local  transmitted  intensity  through  the  negative 

Mn 

and  Iin  being  the  local  incident  intensity  from  the  densitometer. 
A  typical  H&D  curve  is  shown  in  figure  10.  The  interesting  thing  about  that  curve  is  the 
large  linear  region.  If  the  film  is  used  in  this  region  the  photographic  density  can  be  written 
as 

D=ylogE-D0  (16) 

where  y  is  the  slope  of  the  curve  in  the  linear  region  and 
Do  is  the  intersection  of  the  linear  approximation  with  the  D  -  axis. 

If  now  the  definition  of  the  density  D  from  equation  (15)  is  inserted  into  equation  (16),  we 
are  able  to  find  an  expression  for  the  exposure  and  hence  the  intensity  J  in  terms  of  know 

parameters. 

logE  =  -  —  (logx  -  D0) 
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E  =  J  T  =  rO/fl  10<dcVY)  (17) 

The  exposure  time  T  was  measured  for  each  picture,  y  and  Do  were  obtained  through 
some  calibration  negatives.  On  each  roll  of  film  a  couple  of  frames  were  exposed  by 
uniformly  spread  laser  light  at  various  known  exposure  times.  After  scanning  these 
negatives  with  the  densitometer,  each  of  the  negatives  represents  one  point  on  the  actual 
H&D  curve.  Through  these  points  a  linear  least  square  fit  was  made  to  get  y  and  Do.  This 
process  was  done  for  every  roll  of  film,  since  y  and  Do  depend  upon  the  developing 
process. 

After  converting  the  data  negatives  from  relative  densitometer  intensities  to  real 
exposures  and  intensities,  the  angle  calibration  was  applied  to  the  data.  Now  a  first  look  at 
the  data  could  be  taken  (figure  1  la).  In  order  to  check  the  quality  of  our  results  we  have  to 
compare  them  with  valid  models  and  theoretical  calculations.  Therefore  we  go  back  to  Mie 
theory  and  plot  a  Mie  calculation  results  for  each  bubble.  For  the  calculations  the  directly 
measured  radius  of  the  bubble  was  used. 


4.3  Averaging 

By  looking  at  figure  1  la  we  notice  that  the  data  shows  a  lot  more  structure  than  the  Mie 
calculation  result  (figure  12a).  That  suggests  that  the  optical  average  was  not  enough  and 
that  there  is  still  some  noise  in  the  data.  Therefore  we  decided  to  do  running  averages  in 
angle  direction.  First  we  took  an  algorithm  that  averages  the  two  data  points  before  and 
after  with  the  data  point  of  interest;  this  is  called  a  "two  point  average".  Then  we  did  the 
same  but  only  including  one  point  before  and  after  the  data  point  of  interest,  a  "one  point 
average".  The  results  are  shown  in  figures  lib  and  1  lc.  We  notice  two  effects  of  the 
average.  First  of  all  the  data  does  not  look  so  noisy  anymore.  It  now  looks  more  like  the 
Mie  calculation  results  in  figure  12a.  The  second  effect  of  the  averaging  is  that  the 
amplitude  of  the  fringes  is  reduced.  This  would  not  be  bad  because  we  could  compensate 
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Figure  11a:  Intensity  profile  of  experimental  raw  data  obtained  through  a 

microdensitometer  scan.  (A  scan  of  the  same  bubble  shown  in  figure  9, 

ka  =  1278,  a  =  78.5  |im.) 


Figure  lib:  Experimental  data  as  in  figure  1  la,  but  two  point  averaged. 

Figure  11c:  Experimental  data  as  in  figure  1  la,  but  one  point  averaged. 

Figure  12a:  Mie  theory  result  for  the  same  bubble  displayed  in  figure  1  la, 
ka  =  1278,  a  =  78.5  |im,  non  averaged. 

Figure  12b:  Mie  theory  result  as  in  figure  1 2a,  but  two  point  averaged. 

Figure  12c:  Mie  theory  result  as  in  figure  12a,  but  one  point  averaged. 


0.000 


for  that  with  a  scaling  factor  applied  to  the  averaged  data.  But  since  the  amplitude  of  the 
high  frequency  oscillation  is  much  more  affected  by  the  averaging  than  the  low  frequency 
oscillations,  we  get  a  non  uniform  reduction  of  the  amplitude.  This  makes  it  more  difficult 
to  compare  experimental  and  theoretical  results.  Because  of  that  we  took  the  same  averages 
of  the  Mie  theory  data  (figure  12b  and  12c).  The  difference  between  the  two  averaging 
procedures  is  the  following.  The  two  point  average  is  more  effective  in  reducing  the  noise, 
but  it  also  reduces  the  high  frequency  amplitudes  much  more  than  the  one  point  average. 
Because  it  is  most  important  to  analyze  the  structure  of  the  data,  we  decided  to  go  with  the 
two  point  average  for  all  other  data  shown  in  the  appendix.  The  stepsize  for  computing  Mie 
theory  and  averaged  Mie  theory  was  set  to  be  the  same  as  for  the  experimental  data.  Now 
we  are  able  to  compare  experiment  and  theory.  Doing  so,  we  notice  several  things. 

i)  Theory  and  experiment  agree  in  one  very  important  feature,  fringes  caused  by 
interference. 

ii)  The  visibility  of  the  fringe  pattern  is  drastically  reduced  in  the  predicted  region. 

iii)  The  baseline  structure  of  experimental  data  and  Mie  theory  is  different. 

iv)  The  intensity  scale  of  experiment  and  theory  is  different. 

In  the  following  section  each  of  this  points  will  be  discussed  and  interpreted  separately. 
Thereby  we  will  follow  one  specific  bubble  through  the  whole  process  of  analysis.  The 
results  for  more  cases  are  shown  in  the  appendix. 


39 


Chapter  5 


By  comparing  experimental  data  and  theoretical  results  we  see  a  huge  difference 
between  these  two  in  their  overall  shape.  While  the  baseline  structure  of  the  theoretical  data 
appears  to  come  to  a  broad  local  minimum  around  106°,  the  shape  of  the  experimental  result 
is  nothing  like  it.  The  reason  for  that  is  the  following.  The  recorded  intensity  on  the 
negative  In  is  a  composite  of  real  scattered  intensity  from  the  bubble  1^  and  background 


In  =  B  +  Idata  (18) 

Therefore  we  have  to  determine  the  background  B  in  order  to  interpret  the  data.  If  we  take  a 
look  at  the  background  (figure  13),  we  notice  that  the  background  is  a  function  of  scattering 
angle.  Since  the  scattering  event  takes  place  in  water,  a  guess  for  one  background 
component  was  dipole  radiation  of  water. 

When  an  electromagnetic  wave  goes  through  water,  it  forces  the  electrons  of  the  water 
molecules  to  oscillate  with  the  frequency  of  the  light.  Each  water  molecule  then  acts  like  a 
dipole,  radiating  in  all  directions  (figure  14).  Since  the  water  molecules  are  randomly 
positioned,  we  get  scattered  light  from  the  water.  This  phenomena  is  also  known  as 
Rayleigh  scattering.  The  irradiated  intensity  Ip  is  a  function  of  angle  and  can  be  written 


as^.iS 


Id  =  A  cos29 


By  looking  at  the  geometry  of  the  problem  we  observe  that  we  are  just  in  the  position  to  see 
some  dipole  scattering  of  water,  vary  ing  with  angle.  If  we  try  to  fit  the  function  given 
through  equation  19  to  the  background  with  A  as  a  free  strength  parameter  proportional  to 


Figure  14:  Polar  plot  of  In  =  cos2# 

The  scattered  intensity  In  is  displayed  by  the  magnitude  of  the  vector  pointing 
from  the  origin  to  the  arc  starting  at  (1,0)  for  0°  scattering  angle. 


the  exposure  time,  we  find  that  only  the  middle  part  of  the  background  fits  well  (figure  15). 
The  other  parts  are  different  and  because  of  that,  it  is  not  possible  to  use  Rayleigh  scattering 
from  water  to  describe  the  detailed  shape  of  the  background  function. 

Therefore  a  different  approach  was  taken.  Instead  of  subtracting  out  the  Rayleigh 
scattering  of  water,  a  fourth  order  polynomial  was  fitted  to  the  background  data  in  the  range 
of  the  used  data,  e.g.  99°  -  1 12°  (figure  16). 

f(0)  =  a  +  b0  +  c02  +  d03  +  e04  (20) 

This  procedure  provides  a  better  fit  in  the  interesting  region  than  fitting  a  function  to  the 
entire  background  data.  It  is  valid  to  take  only  part  of  the  background  as  long  as  it  includes 
the  range  of  angle  seen  in  the  data  pictures. 

After  determining  a  background  function  this  background  function  f(0)  was  subtracted 
from  the  the  recorded  data  In  to  obtain  the  actual  data  (figure  17). 

Wa  =  In  -  <xf(0)  (21) 

where  a  was  an  adjustable  free  parameter,  a  was  chosen  such  that  the  data  showed  about 
the  same  baseline  structure  as  the  Mie  calculation  results.  One  remaining  problem  however 
is  intensity  matching. 

5.2  Intensity  matching 

The  next  step  is  to  adopt  the  intensities  of  the  background  subtracted  data  to  the 
normalized  intensities  of  the  Mie  theory  results.  All  intensities  shown  so  far  in  figures  of 
raw  data  and  background  subtracted  data  are  basically  in  arbitrary  units.  However  this  does 
not  diminish  the  quality  of  the  data  or  its  validity,  since  the  important  thing  is  the  structure. 
We  now  want  to  make  some  comments  about  problems  and  possible  errors  in  trying  to  find 
normalized  intensities. 

As  described  earlier  the  absolute  intensities  are  obtained  through  the  H&D  curve.  Using 


the  H&D  curve  involves  several  premises. 


The  function  In  =  A  cos:0  fitted  to  the  background. 
Solid  line:  In;  dotted  line:  background 


Fourth  order  polynomial  fitted  to  the  background  between  99°  and  1 12 
Solid  line:  polynomial  fit;  dotted  line:  background 


Figure  16 


Background  subtracted  data  as  in  figure  17a,  but  one  point  averaged. 
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i)  In  producing  the  H&D  curve  uniformly  spread  laser  light  was  used.  To  determine  an  ■ 

absolute  H&D  curve,  the  absolute  incident  intensity  onto  the  negative  would  have  to  be 

known  (equation  14).  However  in  this  experiment  this  could  not  be  directly  measured. 

Therefore  there  is  an  undetermined  scaling  factor.  This  does  affect  the  values  for  y  and  Do. 

ii)  It  was  assumed  that  the  film  was  used  in  the  linear  region  of  the  H&D  curve.  If  | 

however  the  incident  intensity  is  very  low  (or  very  high)  the  actual  H&D  curve  diverges  J 

from  linear  behavior  (figure  10).  Since  the  linear  approach  is  used  throughout  all  data,  this  ^ 

i 

might  be  a  source  of  error  especially  for  low  intensities.  Inspection  of  figure  10  suggests  j 

that  one  effect  of  nonlinearity  of  the  H  &  D  curve  would  be  to  offset  the  inferred  intensities. 

In  our  final  analysis  we  add  a  constant  to  the  intensities  in  order  to  compensate  for  these 

and  other  effects.  The  value  for  this  constant  is  taken  to  be  a  free  parameter  for  each .  ! 

! 

negative. 

iii)  Furthermore  it  was  assumed  that  the  densitometer  reading  is  a  linear  function  of  the 

incident  light  provided  by  that  device.  However  this  may  not  be  exactly  the  case.  The  ! 

i 

output  of  the  densitometer  may  be  slightly  nonlinear.  However  this  behavior  does  not  affect 
the  general  shape  and  structure. 

iv)  The  air  glass  water  interfaces  at  the  entrance  window  and  observation  window  were 
not  accounted  for.  Although  they  are  perpendicular  to  most  of  the  light,  they  do  introduce 
an  error  in  obtaining  intensities  which  might  depend  weakly  on  angle.  The  severity  of  the 
problem  is  reduced  because  of  the  antireflection  coating  on  the  outside  of  the  observation 
window.  The  coating  reduces  the  strength  of  internal  reflections  slightly.  Additional 
experiments  on  scattering  from  a  reflecting  cylindrical  wire  suggest  that  the  transmittance 
through  the  observation  window  depends  very  weakly  on  angle. 

Therefore,  the  following  procedure  to  match  experimental  and  theoretical  intensities  was 
applied.  Since  all  theoretical  results  are  normalized  to  the  iiradiance  of  a  perfectly  reflecting 
sphere  of  the  same  size  as  the  bubble,  the  experimental  data  were  multiplied  by  a  scaling 
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factor  (3  such  that  the  distance  between  minima  and  maxima  of  fringes  at  a  given  angle  were 
the  same  as  for  the  Mie  theory  result  of  that  particular  bubble.  Then  a  constant  c  was 
subtracted  to  lower  the  whole  spectrum  to  the  desired  level  defined  by  the  Brewster  angle 
scattering  region.  Both  constants,  for  multiplication  and  subtraction,  were  necessary 
because  of  the  above  described  implications  in  obtaining  absolute  intensities. 

^Norm  =  P  Idata  "  C  (22) 

The  final  result  is  shown  in  figure  18.  These  results  are  now  taken  into  further 
investigations  about  fringe  patterns  and  their  agreement  with  the  theory. 


5.3  Disappearing  of  fringe  pattern 

The  clearly  visible  reduction  of  the  fringe  pattern  strongly  suggests  that  we  do  see  the 
Brewster  angle  effect  in  the  experiment.  We  especially  note  the  disappearing  of  the  high 
frequency  oscillation  which  was  attributed  to  interference  of  the  (0,0)  and  (2,1)  rays.  That 
indicates  that  the  (0,0)  ray  really  goes  through  a  local  minimum  in  intensity.  The  (2,1)  ray 
however  is  still  there  and  causes  low  angular  frequency  interference  with  the  (3,2)  ray  as 
identified  earlier.  We  also  find  the  observed  location  of  the  Brewster  angle  region  in  general 
agreement  with  the  calculated  results  from  both,  equation  6  and  the  Mie  theory  data  shown 
in  figure  12.  By  comparing  figures  12  and  18  we  note  that  the  general  structure  in  the 
Brewster  angle  region  is  quite  similar.  This  is  a  very  nice  result,  considering  the  weakness 
of  the  signal. 

However  there  is  one  source  of  error.  We  always  tried  to  launch  the  bubbles  such  that 
they  rose  through  the  center  of  the  laser  beam.  However  due  to  the  statistical  nature  of  that 
process,  it  was  not  always  possible  for  the  bubble  to  go  straight  through  the  center  of  the 
beam.  The  problem  in  that  case  is  that  the  relative  magnitudes  of  the  rays  contributing  to  the 
scattering  pattern  could  be  out  of  proportion.  Because  all  relevant  rays  enter  the  bubble  at  a 
different  point  those  rays  closer  to  the  edge  of  the  beam  are  weaker  than  they  should  be  and 


-  '  v  v ' 


18a:  Final  result:  background  subtracted,  two  point  averaged,  scaled  data, 

ka  =  1278,  a  =  78.5  pm.  The  85  mm  lens  was  used.  The  scaling  was  done 
according  to  the  non  averaged  Mie  theory  result  of  figure  12a. 

18b:  Final  result,  the  same  experimental  data  as  in  figure  1 8a,  but  scaled  to  the  two 
point  averaged  Mie  theory  result  of  figure  12b. 


18c:  Final  result,  the  same  experimental  data  as  in  figure  1 8a,  but  scaled  to  the  one 
point  averaged  Mie  theory  result  of  figure  12c. 


the  rays  closer  to  the  center  of  the  beam  are  relatively  strong.  Therefore  we  may  not  always 
find  the  relative  magnitude  of  the  fringes  to  be  as  in  the  Mie  theory  results.  Fortunately  the 
bubble  diameter  was  always  much  smaller  than  the  beam  width,  so  that  this  effect  should 
usually  be  small. 

5.4  Fringe  pattern 

Both  experiment  and  theory  exhibit  the  predicted  oscillatory  fine  structure.  As  shown 
earlier  this  structure  is  caused  by  interference  from  mainly  the  (0,0)  ray  with  the  (2,1)  ray. 
To  check  results  of  the  experiment  we  calculate  the  separation  between  adjacent  extreme  A0 
for  a  given  measured  size  using  formula  11.  Comparing  fringe  spacing  between  theoretical 
and  experimental  results  gives  information  about  the  size  of  the  bubble  since  the  separation 
between  adjacent  extreme  is  inversely  proportional  to  the  size.  Results  of  this  comparison 
are  shown  in  table  4.  The  numbers  were  obtained  by  superposing  a  very  fine  grid  on  the 
plot  and  using  it  to  measure  the  period  of  interest.  Using  the  measured  bubble  radius,  the 
fringe  spacing  for  all  possible  interferences  was  calculated.  Then  a  comparison  with  the 
fringe  spacing  detected  in  the  experimental  data  was  done.  In  addition  from  the  observed 
fringe  spacing,  a  radius  value  for  each  bubble  was  calculated  using  the  measured  (0,0), 
(2,1)  interference  period.  This  is  the  most  pronounced  structure  and  in  all  cases  the 
(0,0),(2,1)  interference  could  be  measured  with  good  accuracy.  This  inferred  radius  value 
was  then  compared  with  the  measured  radius.  It  was  found  to  be  in  general  agreement 
considering  the  circumstances  of  the  size  measurement.  The  measuring  microscope  had  an 
intrinsic  error  of  about  ±3  Jim.  Besides  that  there  was  always  the  problem  of  determining 
the  edge  of  the  bubble  while  trapped  under  the  glass  slide  and  measured.  Furthermore  there 
we  had  to  deal  with  the  problem  of  vibration.  Considering  the  size  of  the  bubbles,  even 
very  small  vibrations  could  add  a  considerable  error  to  the  measurement.  All  together  we 
believe  to  have  an  error  of  about  ±10  pm  in  the  measurement  of  the  bubble  diameter. 
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Table  4:  Analysis  of  the  occurring  interferences 


Radius  of  the  bubble  in  microns 
observed  calculated 

78.5 0000  76.3 


Interference  Angular  separation  in  deg. 


calculated 


observed 


Difference 


-2.2 


Difference 


(0,0)  (2,1) 


0.291768 


2 .7% 


(0,0)  (3,1) 

(0,0)  (3,2) 


0.232058 


0.936399 


0.85 


9.2% 


(2,1)  (3,1) 

(2,1)  (3,2) 


1.133936 


0.423825 


(3,1)  (3,2) 


0.308514 


Inspection  of  the  data  reveals  a  low  frequency  modulation  outside  the  Brewster 
scattering  angle  region  caused  by  (0,0), (3, 2)  interference.  Right  at  the  Brewster  scattering 
angle  region,  basically  only  the  (2,1), (3, 2)  interference  was  found.  However  a  problem  in 
this  region  is  the  weakness  of  the  signal.  Therefore  the  signal  is  much  more  succeptible  to 
noise,  which  makes  it  a  lot  harder  to  determine  the  frequency  of  the  fringe  spacing.  Another 
source  of  error  is  the  earlier  discussed  size  measurement  of  the  bubble.  Since  the  fringe 
spacing  is  very  sensitive  to  the  radius  of  the  bubble,  inaccuracies  in  determining  the  size  of 
the  bubble  cause  variations  in  A8.  Nevertheless  we  get  fairly  good  agreement  between 
theory  and  experiment. 
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Chapter  6 


The  principal  effects  of  Brewster  angle  reflection  on  the  scattering  from  uncoated  air 
bubbles  in  water  were  observed  experimentally.  The  final  results  also  show  all  the 
important  features  predicted  by  Mie  calculations.  The  disappearing  fringe  pattern  as  well  as 
the  relative  minimum  in  scattered  intensity  were  observed  experimentally.  The  periods  of 
interference  patterns  for  several  refracted  and  reflected  rays  were  found  to  agree  reasonably 
well  with  the  measured  size  of  the  bubble  and  geometrical  consideration.  These  effects 
might  be  used  to  distinguish  between  bubbles  and  other  scatterers7-9.  (Solid  particles  in 
water  would  not  show  a  Brewster  angle  effect  near  106°.)  Furthermore  this  technique  might 
be  useful  to  detected  coatings  on  bubbles  since  the  Brewster  angle  effect  is  strongly 
dependent  upon  the  relative  index  of  refraction.  If  investigated  more,  it  might  become  a 
useful  tool  for  getting  information  about  the  type  of  coating  material  and  the  thickness  of 
the  coat  on  bubbles  in  dielectric  media7-9. 
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APPENDIX  A 


In  the  Appendix  a  few  more  examples  of  experimental  data  shall  be  presented.  For  each 
case  the  experimental  data  is  shown  next  to  the  Mie  theory  results.  This  is  done  with  and 
without  the  two-point  average.  Furthermore  a  complete  analysis  of  all  detectable  fringes 
was  done.  The  results  were  obtained  in  the  same  way  as  described  earlier. 


Table  Al:  Analysis  of  the  occurring  interferences 


Radius  of  the  bubble  in  microns 


observed 

61.00000 


calculated 

65.4 


Interference  Angular  separation  in  deg. 


calculated 


observed 


Difference 


+  4 . 4 


Difference 


(0,0)  (2,1) 

(0,0)  (3,1) 

(0,0)  (3,2) 
(2,1)  (3,1) 

(2.1)  (3,2) 

(3.1)  (3,2) 


0.375472 


0.298632 


1.205037 


1.459246 


0.545415 


0.397022 


0.35 


0.52 


6.8% 


8.7% 


4.7% 


Radius  of  the  bubble  in  microns 


observed 

80.00000 


calculated 

76.3 


Difference 


Interference  Angular  separation  in  deg. 


calculated 


observed 


Difference 


(0,0)  (2,1) 


0 . 286297 


4.6% 


(0,0)  (3,1) 


0.227707 


(0,0)  (3,2) 


0.918841 


2  .  1% 


(2,1)  (3,1) 

(2,1)  (3,2) 


1.  112675 


0 . 415879 


(3,1)  (3,2) 


0.302729 
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Figure _ A3 

Radius  of  the  bubble  in  microns  Difference 

observed  calculated 


72 

.  00000 

76.3 

+  4.3 

Interference 

Angular  separation 
calculated 

in  deg. 
observed 

Difference 

(0,0) 

(2,1) 

0.318108 

0.3 

5.7% 

(0,0) 

(3,1) 

0.253008 

(0,0) 

(3,2) 

1.020935 

0.85 

16.7% 

(2,1) 

(3,1) 

1.236306 

(2,1) 

(3,2) 

0.462087 

0.45 

2.6% 

(3,1) 

(3,2) 

0.336366 

Ei-gucs _ M 

Radius  of  the 
observed 
73.50000 

bubble  in  microns 
calculated 

65.4 

Difference 

-8  .  1 

Interference 

Angular  separation 
calculated 

in  deg. 
observed 

Difference 

(0,0) 

(2,1) 

0.311616 

0.35 

11.0% 

(0,0) 

(3,1) 

0.247844 

(0,0) 

(3,2) 

1.000099 

1.025 

2 . 5% 

(2,1) 

(3,1) 

1.211075 

(2,1) 

(3,2) 

0. 452657 

0 . 55 

17.3% 

(3,1)  (3,2) 


0. 329501 
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Figure  A5 


Radius  of  the 
observed 

71 . 50000 

bubble  in  microns 
calculated 

83 . 4 

Interference 

Angular  separation 
calculated 

in  deg. 
observed 

(0,0)  (2,1) 

0.320332 

0.275 

(0,0)  (3,1) 

0.254777 

(0,0)  (3,2) 

1.028074 

0.825 

(2,1)  (3,1) 

1.244951 

(2,1)  (3,2) 

0.465319 

0.45 

(3,1)  (3,2) 

0.338718 

Fiqure  A6 

Radius  of  the 
observed 

61 . 50000 

bubble  in  microns 
calculated 

65.4 

Interference 

Angular  separation 
calculated 

in  deg. 
observed 

(0,0)  (2,1) 

0.372419 

0. 35 

(0,0)  (3,1) 

0.296204 

(0,0)  (3,2) 

1. 195240 

1.1 

(2,1)  (3,1) 

1.447382 

(2,1)  (3,2) 

0. 540980 

0.55 

(3,1)  (3,2) 

0 . 393794 

+  11.9 

Difference 

14.2% 

19.8% 

3 .3% 

Difference 
+  3.9 

Difference 

6.0% 

8 . 0% 

1 . 6% 


Experimental  data,  ka  =  993,  a  =  61  jim,  135  mm  lens,  two  point  averaged 
and  scaled  to  the  non  averaged  Mie  theory  result. 

Mie  theory  result,  ka  =  993,  non  averaged 

Experimental  data  as  in  figure  A  1.1  a,  but  scaled  to  the  averaged  Mie  theory 
result. 

Mie  theory  result  as  in  figure  Al.2a,  but  two  point  averaged. 
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L  Experimental  data,  ka  =  1303,  a  =  80  (im,  135  mm  lens,  two  point 
averaged,  scaled  to  the  non  averaged  Mie  theory  result 


i2.2a:  Mie  theory  result,  ka  =  1303,  non  averaged 


2i  Experimental  data  as  in  figure  A2.1a,  but  scaled  to  the  averaged  Mie  theory 
result. 


>.2b:  Mie  theory  result  as  in  figure  A2.2a,  but  two  point  averaged. 


Figure  A2.1b 
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Experimental  data,  ka  =  1 172,  a  =  72  (am,  85  mm  lens,  two  point  averaged, 
scaled  to  the  non  averaged  Mie  theory. 

Mie  theory  result,  ka  =  1172,  non  averaged 

Experimental  data  as  in  figure  A3. la,  but  scaled  to  the  averaged  Mie  theory 
result. 

Mie  theory  result  as  in  figure  A3. 2a,  but  two  point  averaged. 


Figure  A3. la 
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Experimental  data,  ka  =  1 197,  a  =  73.5  (im,  85  mm  lens,  two  point 
averaged,  scaled  to  the  non  averaged  Mie  theory  result. 

Mie  theory  result,  ka  =  1 197,  non  averaged. 


Experimental  data  as  in  figure  A4.1a,  but  scaled  to  the  averaged  Mie  theory 
result. 


Mie  theory  result  as  in  figure  A4.2a,  but  two  point  averaged. 


Figure  A4. lb 
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Figure  A4.2s 
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Experimental  data,  ka  =  1 164,  a  =  71.5  (im,  85  mm  lens,  two  point 
averaged,  scaled  to  the  non  averaged  Mie  theory. 

Mie  theory  result,  ka  =  1164,  non  averaged. 

Experimental  data  as  in  figure  A5. 1  a,  but  scaled  to  the  averaged  Mie  theory 
result. 

Mie  theory  result  as  in  figure  A5.2a,  but  two  point  averaged. 
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